Stellar and gaseous abundances in M82 
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ABSTRACT 

The near infrared (IR) absorption spectra of starburst galaxies show several 
atomic and molecular lines from red supergiants which can be used to infer reli- 
able stellar abundances. The metals locked in stars give a picture of the galaxy 
metallicity prior to the last burst of star formation. The enrichment of the new 
generation of stars born in the last burst can be traced by measuring the hot 
gas in the X-rays. For the first time detailed stellar abundances in the nuclear 
region of the starburst galaxy M82 have been obtained. They are compared with 
those of the hot gas as derived from an accurate re- analysis of the XMM- Newton 
and Chandra nuclear X-ray spectra. The cool stars and the hot gas suggest 
[Fe/H]= —0.35 ± 0.2 dex, and an overall [Si,Mg/Fe] enhancement by ~0.4 and 
0.5 dex, respectively. This is consistent with a major chemical enrichment by 



^Based on observations with the Near Infrared Camera Spectrometer (NICS) mounted at the Telescopic 
Nazionale Gahleo (TNG), La Palma, Spain and archival XMM- Newton and Chandra data. 
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SNe II explosions in recursive bursts on short timescales. Oxygen is more puz- 
zling to interpret since it is enhanced by ^0.3 dex in stars and depleted by 
~0.2 dex in the hot gas. None of the standard enrichment scenarios can fully 
explain such a behavior when compared with the other ct-elements. 

Subject headings: galaxies: individual (M82) — galaxies: abundances — galaxies: 
starburst — X-ray: galaxies — infrared: galaxies 

1. Introduction 

The signature of the star formation (SF) history of a galaxy is imprinted in the abun- 
dance patterns of its stars and gas. Determining the abundance of key elements released in 
the interstellar medium (ISM) by stars with different mass progenitors and hence on different 
time scales, will thus have a strong astrophysical impact in drawing the global picture of 
galaxy formation and evolution (McWiUiam 1997). 

So far, the metallicities of distant starburst (SB) galaxies have been mainly derived by 
measuring the nebular lines associated to their giant Hll regions (Storchi-Bergmann, Calzetti 
& Kinney 1994; Coziol et al. 1999; Henry, Edmunds & Koppen 2000). Direct determination 
of O and N abundances can be obtained once electron temperature (Tg) is well estabhshed 
(McCall 1984). At high metallicities Tg decreases due to strong lines cooling and its major 
diagnostic (the weak [OIII] A4636) is no longer observable. Several alternative methods based 
on other strong lines were proposed in the past years (Pagel et al. 1979; Dopita & Evans 1986; 
McGaugh 1991) to derive Tg, but all of them rely on a large number of physical parameters 
and model assumptions for the geometry and the nature of the ionizing populations (c.f. 
c.g Stasinska (2001) for a review). Depending on the adopted calibrations, quite different 
abundance sets can reproduce the observed line ratios. 

However, abundances in SB galaxies can be also measured from absorption stellar fea- 
tures in the near IR and/or emission lines in the hot X-ray gas. Stellar abundances can 
be poorly constrained from optical spectra, since the nebular emission strongly dilutes the 
absorption hues and dust can heavily obscure the central regions where most of the burst 
activity is concentrated. On the contrary, in the near IR the stellar continuum due to red 
supergiants (RSGs) usually largely dominates over the possible gas and dust emission (Oliva 
& Origha 1998; Origha & Oliva 2000). The X-ray emission from SB galaxies at E>2 keV 
is mainly due to high mass X-ray binaries with power-law spectra, while at lower energies 
hot plasma mainly heated by SN explosions (Dahlem et al. 1998) dominates and the re- 
sulting spectrum is a superposition of thermal continuum and emission lines from atomic 
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recombinations (Fabbiano 1989; Persic & Rephaeli 2002). 

Determining the abundances of SB galaxies is not only relevant to constrain their SF 
history, but also offers the unique chance of directly witnessing the enrichment of the ISM 
(see e.g. Maeder & Conti 1994, for a review). Metals locked into stars give a picture of the 
enrichment just prior to the last burst of SF, while the hot gas heated by SNe II explosions 
should trace the enrichment by the new generation of stars. Nebular metallicities from the 
cold gas should potentially trace the stellar ones if the cooling and mixing timescales are 
slow (up to 1 Gyr), as suggested by some models (see e.g. Tenorio-Tagle 1996) since the 
gas enriched by the ongoing burst is still too hot to be detected in the optical lines. On 
the contrary, if rapid (a few Myr) cooling and mixing occur as suggested by other models 
(see e.g. Recchi, Matteucci & D'Ercole 2001), the nebular abundances could trace the gas 
enriched by the new population of stars (similarly to X-rays). Hence, important constraints 
on the cooling and mixing timescales of the gas can be provided by comparing the metallicity 
inferred from stars and hot plasma with those of the cold gas. 

Sect. 2 briefly reviews the nebular abundance estimates of M82, as published so far in 
the literature. Sect. 3 shows our near IR spectra of the nuclear region of M82 and the derived 
stellar abundances. Sect. 4 presents a re-analysis of the XMM-A^eu;ton and Chandra nuclear 
X-ray spectra and the hot gas-phase abundances. Sect. 5 summarized the overall chemical 
enrichment scenario in M82 as traced by the stellar, cold and hot gas-phase abundances, 
while in Sect. 6 we draw our conclusions. 



2. Nebular abundances of M82 from the literature 

M82 is a prototype of SB galaxies (Rieke et al. 1980), experiencing a major SF formation 
episode in its nuclear regions, with strong super- wind and SN activity (Wills et al. 1999). 
Despite the many multi-wavelength studies (see e.g. Golla, Allen & Kronberg 1996, and 
references therein), only few measurements of nebular abundances, mainly derived from 
mid/far IR forbidden lines, have been published so far. From the analysis of the [Sill] 35fim 
line in the central 300 pc in radius at the distance of M82 (i.e. d~3.6 Mpc), Lord et al. 
(2001) suggest an indicative [Si/H] ~ —0.4 dex relative to the Solar value of Grevesse & 
Sauval (1998). 

More recently, Forster Schreiber et al. (2001), by analyzing the nuclear spectra of M82 
taken with the Short Wavelength Spectrometer (SWS) onboard the Infrared Space Observa- 
tory (ISO) within an aperture corresponding to 100 x 200 pc on the sky, give cold gas-phase 
abundances of [Ne/H]=-F0.1, [Ar/H]=-F0.3 and [S/H]=-0.7 dex relative to the Solar values 
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of Grevesse & Sauval (1998), with an overall uncertainty between 0.2 and 0.3 dex. 

By taking the hne ratios listed by AUoin, CoUin-Souffrin, & Joly (1979) which refer 
to the central ~100 pc in radius of M82 (Peimbert & Spinrad 1970), and the empirical 
cahbrations of McGaugh (1991) which give N and Fe abundances as a function of [0/H] 
(see e.g. their Figure 11, equations 9 and 10a), one can also obtain rough estimates of 
[O/H]=+0.1, [N/H]=0.0 and [Fe/H]=-0.3 dex relative to the Solar values of Grevesse & 
Sauval (1998), with an overall uncertainty of ±0.3 dex. 

The first attempt to measure hot gas-phase abundances in M82 has been made quite 
recently by Ptak et al. (1997) and slightly revised by Umeda et al. (2002), by analyzing ASCA 
X-ray (0.4-10.0 kcV) spectra. Fc and O abundances as low as 1/20 Solar and significantly 
higher (by a factor between 3 and 10) Si, S, Mg, Ca ones have been obtained, also consistent 
with BeppoSAX results (Cappi et al. 1999). 

More recently. Read & Stevens (2002) by analyzing the spectra obtained with the Re- 
flection Grating Spectrometer (RGS) onboard XMM-A'^etf^ton found near-Solar Fe and O and 
super-Solar (between a factor of 2 and 5) Mg, Si, Ne, N abundances. While the abundance 
ratios are in reasonable agreement within ±0.3 dex with Ptak et al. (1997), there is one order 
magnitude difference between the zero-point of the two calibrations. 

3. IR spectra and stellar abundances 

Near-IR spectroscopy is a fundamental tool to obtain accurate abundances of key 
elements like Fe, C, O and other metals (e.g. Si, Mg, Ca, Na, Al, Ti) in cool stars 
(Teff < 5000 K). Several atomic and molecular lines arc strong, not affected by severe 
blending, hence also measurable at low- medium resolution (Kleinmann & Hall 1986; Origlia 
et al. 1993; Wallace & Hinkle 1997; Joyce et al. 1998; Frogel et al. 2001), making them 
powerful abundance tracers not only in stars but also in more distant stellar clusters and 
galaxies and in a wide range of metaUicities and ages (Origha et al. 1997; Ohva & Origlia 
1998). Abundance analysis from integrated spectra of galaxies requires full spectral synthe- 
sis techniques to properly account for line blending and population synthesis to define the 
dominant contribution to the stellar luminosity. 

In SB galaxies the stellar IR continuum usually dominates over the possible gas and 
dust emission. This represents a major, conceptual simplification in population and spectral 
synthesis techniques, making possible and easier the interpretation of integrated spectra from 
distant stellar clusters and galaxies. 




Fig. 1. — Near-IR spectra of the nuclear region of M82. Observed spectra: histograms; 
synthetic stellar best-fit solution and two other models with ±0.3 dcx abundance variations: 
solid lines. A few major stellar and nebular lines are also marked. 
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Since the near IR continuum of SB galaxies is dominated by luminous RSGs (see e.g. 
Origlia & Oliva 2000, for a review), their integrated spectrum can be modeled with an 
equivalent, average star, whose stellar parameters (temperature Teff, gravity logg and mi- 
croturbulence velocity ^) mainly depend on the stellar age and metallicity. Both observations 
and evolutionary models (see e.g. Keller 1999; Origlia et al. 1999; Origlia 2003, and reference 
therein) suggest that RSGs of ages between ~6 and 100 Myr and metallicities between 1/10 
and Solar are characterized by low gravities (log g<1.0), low temperatures (< 4000 K) and 
relatively high microturbulence velocity (^ >3 km/s). A variation in the adopted stellar pa- 
rameters for the average RSG population by ATeff±200 K, Alog g=±0.5 and A^^O.5 km 
implies a < ±0.1 dex change in the abundances estimated from atomic lines, and < ±0.2 dex 
in those estimated from the molecular lines, which are more sensitive to stellar parameters. 

1.0-1.8/i longslit spectra (see Fig. 1) at R=2,500 of the nuclear region of M82 were 
obtained with the Near IR Camera Spectrometer (NIGS) mounted at the Itahan Telescopio 
Nazionale Galileo (TNG) on December 2002. The spectra were sky-subtracted, fiat-fielded 
and corrected for atmospheric absorption using an 0-star spectrum as reference. They 
were wavelength calibrated by using a Ne-Ar lamp and the mono dimensional spectra were 
extracted by summing over the central 0.5" x 3.0", corresponding to an aperture projected 
on the source of 9 x 52 pc at the distance of M82. The spectra have been normalized to the 
continuum, which was determined applying a low-pass smoothing filter to each spectrum. 
Total integration times of 72 and 32 min in the J and H bands, respectively, were used, 
providing a final signal to noise ratio >40. 

By measuring the absorption line broadening a stellar velocity dispersion a ~ 105±20 
km/s has been derived, in perfect agreement with the estimate by Gafney, Lester & Telesco 
(1993) and in the typical range of values measured in other massive SB galaxies (Ohva et 
al. 1999). 

A grid of synthetic spectra of red supergiant stars for different input atmospheric 
parameters and abundances have been computed, using an updated (Origlia et al. 2002, 
2003) version of the code described in Origlia et al. (1993). Briefly, the code uses the 
LTE approximation and is based on the molecular blanketed model atmospheres of John- 
son, Bernat & Krupp (1980). It includes several thousands of near IR atomic lines and 
molecular roto- vibrational transitions due to CO, OH and CN. Three main compilations 
of atomic oscillator strengths are used, namely the Kurucz's database (c.f. http://cfa- 
www.harward.edu/amdata/ampdata/kurucz23/sekur. html), and those published by Biemont 
& Grevesse (1973) and Melendez k Barbuy (1999). 

Our code provides full spectral synthesis over the 1-2.5 /um range. Given the high 
degree of line blending, our abundance estimates are mainly obtained by best-fitting the full 
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Fig. 2. — Left panel: probability of a random realization of our best-fit full spectral synthesis 
solution with varying the elemental abundances A[X/H] of ±0.2 and 0.3 dex with respect to 
the best-fit (see Sect. 3). Right panel: average difference between the model and the observed 
equivalent widths of a few selected hues (see Sect. 3 and Fig. 1). 
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observed spectrum and by measuring the equivalent widths of a few selected features (cf. 
Fig. 1), dominated by a specific chemical element, as a further cross-check. The equivalent 
widths have been measured by performing a Gaussian fit with a equal to the measured stellar 
velocity dispersion, typical values ranging between 0.5 and 3 A, with a conservative error of 
±200 mA to also account for a ±2% uncertainty in the continuum positioning. 

By fitting the full observed IR spectrum and by measuring the equivalent widths of se- 
lected lines, we obtained the following best-fit stellar parameters and abundance patterns for 
M82: Tej/=4000, log5=0.5, ^=3, [Fe/H]=-0.34; [O/Fe]=+0.34, [<Si,Mg,Ca>/Fe]=+0.38; 
[Al/Fe]=+0.5; [C/Fe]=-0.26; i2c/i3c<10. Table 3 lists the derived abundances and their 
associated random errors at 90% confidence. Reference Solar abundances are from Grevesse 
& Sauval (1998). 

Synthetic spectra with lower element abundances are systematically shallower than the 
best-fit solution, while the opposite occurs when higher abundances are adopted. In order to 
check the statistical significance of our best-fit solution, as a function of merit we adopt the 
difference between the model and the observed spectrum (hereafter S). In order to quantify 
systematic discrepancies, this parameter is more powerful than the classical test, which is 
instead equally sensitive to random and systematic scatters (see Origlia et al. 2003, for more 
details) . 

Since 5 is expected to follow a Gaussian distribution, we compute 5 and the corre- 
sponding standard deviation (a) for the best-fit solution and 6 test models with abundance 
variations A[X/H] — ±0.1, 0.2 and 0.3 dex with respect to the best-fit. We then extract 
10000 random subsamples from each test model (assuming a Gaussian distribution) and we 
compute the probability P that a random realization of the data-points around a test model 
display a d that is compatible with the best-fit model. P ~ 1 indicates that the model is a 
good representation of the observed spectrum. 

The left panel of Fig. 2 shows the average results for the observed J and H band spectra 
of M82. It can be easily appreciated that the best-fit solution provides in all clear 
maximum in P (>99%) with respect to the test models. More relevant, test models with an 
abundance variation A[X/H] > ±0.2 dex lie at ~ 1.5a from the best-fit solution, while test 
models with A[X/H] > ±0.3 dex lie at u > 2 from the best-fit solution. 

The analysis of the line equivalent widths provide fully consistent results. The right 
panel of Fig. 2 shows the average difference between the model and the observed equivalent 
width measurements. Models with ±0.2 dex abundance variations from the best-fit solution 
are still acceptable at a ~ 1.5a significance level, while those with ±0.3 dex variations are 
only marginally acceptable at a 2 — 3(7 level. 
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Tablc 1. Element abundances in Solar units" as 
derived from our analysis of the IR and X-ray spectra. 
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Note. — Solar values from Grevesse & Sauval (1998). 
^ Stellar abundances from TNG/NICS near IR absorption 
spectra of red supergiants (see Fig. 1 and Sect. 3). 

Hot gas- phase abundances from XMM/RGS and pn 
spectra (see Figs. 3,5 and Sect. 4). 
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Models with stellar parameters varying by ATeff±200 K, Alog g=±0.5 and A^^O.5 km s~^ 
and abundances varying accordingly by 0.1-0.2 dex, in order to still reproduce the deepness 
of the observed features, are also less statistical significant (on average only at > 2a level) 
with respect to the best-fit solution. Hence, as a conservative estimate of the systematic er- 
ror in the derived best-fit abundances, due to the residual uncertainty in the adopted stellar 
parameters, one can assume a value of < ±0.1 dex. 

By taking into account the overall uncertainty in the definition of the average population 
and the statistical significance of our spectral synthesis procedure, we can safely conclude that 
the stellar abundances can be constrained well within ±0.3 dex and their abundance ratios 
down to ~0.2 dex, since some (if not all) of the stellar parameter degeneracy is removed. 

4. X-ray spectra and hot gas-phase abundances 

The determination of the abundances of the hot X-ray emitting gas in SB galaxies has 
traditionally suffered from large uncertainties. Indeed, the low angular and spectral resolu- 
tion of the various X-ray telescopes in the prc-XMM-A^ett;ton/ Chandra era did not allow to 
disentangle between point sources and hot gas emission, making abundance determinations 
severely model-dependent (Dahlem et al. 2000). The problem hes in the fact that the X-ray 
spectrum of SB galaxies contains two major components: the emission from hot diffuse gas 
and the integrated contribution of point sources. The first is described by an optically thin 
thermal spectrum plus emission lines, the latter by a power-law. If the angular resolution is 
not good enough, it is not possible to reliably subtract the point sources from the total spec- 
trum, so that the equivalent widths of the emission lines (and thus the element abundances) 
are not unambiguously defined. 

On the other hand, high angular resolution alone is not enough: recent Chandra stud- 
ies (Strickland ct al. 2002; Martin, Kobulnicky & Hcckman 2002) demonstrated that the 
relatively low spectral resolution of the ACTS detector makes individual element abundance 
analysis still problematic. The high spectral resolution of the Chandra gratings rapidly 
degrades for sources which are much more extended than the instrumental Point Spread 
Function (PSF) and thus these gratings are almost useless for the study of nearby SB galax- 
ies. However, by coupling the high angular resolution (0.5 PSF) of Chandra with the high 
spectral resolution (A/AA ~ 400) of the XMM/RGS, which is not sensitive to the source 
extension, it is possible to overcome the above described difficulties. 

M82 was observed several times by Chandra; we consider the two longest exposures, of 
33 ks and 15 ks respectively, both obtained on September 1999. A cumulative spectrum of 
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the brightest point sources was extracted from the two Chandra ACIS-I observations; the 
best-fit model was an absorbed power-law with Nh = 7.9 ±0.7 x 10^^ cm~^, F = 0.84 ±0.07, 
Xr = 0.4, which is typical of High Mass X-ray Binaries (see Persic & Rephaeli 2002, and 
reference therein). The 0.5-10 kcV flux was 4.2 x 10^^^ erg s"^ cm~^, which represented 
~ 25% of the total (point source plus diffuse) flux (but notice that M82 is a variable source, 
Ptak & Griffiths 1999, Rephaeli & Gruber 2002). 

M82 was observed by XMM- Newton for 30 ks in May 2001; in the XMM- Newton ob- 
servation, after light-curve cleaning, ^ 20 ks of data were available for scientiflc analysis. 
Wc extracted the EPlC-pn spectrum from a circular region centered on the starbursting core 
of M82 and with a 15" radius, which approximately matches the RGS PSF; the standard 
recipe for EPIC spectra extraction suggested in the XMM web pages was followed. The 
XMMSAS software version 5.3.0 was used. A background spectrum was taken from standard 
background files provided by the XMM-A^eu;ton SOC. RGS spectra were extracted with stan- 
dard settings (i.e. including events within 90% confidence and rejecting as background events 
those outside 95% of the PSF figure); this procedure is similar to the one used by Read & 
Stevens (2002) but with more recent calibrations. After visual inspection, we discarded the 
2nd-order spectra because of low signal/noise ratio. We included in the fit only the Ist-order 
channels comprised in the 6-20 A interval, since outside the S/N ratio rapidly drops. A sys- 
tematic error up to 10% in the 14-18 A channels was included (XMM cahbration document 
XMM-SOC-CAL-TN-0030 Issue 2). 

The hard-band (2-7 keV) Chandra images show that the point sources emission largely 
exceeds that of the hot gas, while the hard part (3-7 keV) of the pn spectrum can be 
well described by a power-law whose best-fit slope resulted in the same value previously 
determined for the Chandra point sources. Thus, making use of information from both 
Chandra and XMM-A^eu;ton observations, we have been able to constrain the point source 
slope, absorption and normalization at the moment of the XMM-A^eu;ton observation. 

We jointly fitted the pn and RGS spectra with XSPEC Version 11.2 (see Fig. 3) with 
a two component model which is the sum of: i) an absorbed power- law (accounting for 
point sources) with the slope and absorption parameters fixed at the best fit values obtained 
with Chandra while normalization was left free to account for variability; ii) an absorbed, 
optically-thin thermal plasma emission, with variable line intensities and differential emission 
measure (DEM) distribution. The DEM distribution was modeled with a 6^^ order Cheby- 
shev polynomial (c6pvmkl model in XSPEC); with the best-fit parameters it resembles a 
bell-shaped curve with a peak at kT ~ 0.7 keV and FWHM ~ 0.65 keV (i.e. the plasma 
temperatures range from ~0.35 keV to ~1 keV), in agreement with Read & Stevens (2002). 
The thermal component is absorbed by cold gas with best-fit column density Nh — 3.8 x 10^^ 



- 12 - 




10 15 



channel wavelength (A) 

Fig. 3. — XMM-RGS spectrum of the nuclear region of M82. Upper panel: data points 
(lower spectrum - open circles: RGSl; filled triangles: RGS2; upper spectrum: pn) along 
with the best-fit model. Lower panel: residuals in units of a. 
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cm ^. 

The confidence contours of the Iron abundance versus normahzation of the thermal 
component are reported in Fig. 4. As expected the two parameters are well correlated as 
the higher is the normalization of the continuum spectrum the lower is the intensity of 
the emission lines and in turn the element abundances. The derived abundances and their 

associated 90% confidence error arc shown in Table 3. We notice that these abundances are 
about a factor 2 below those reported in Read & Stevens (2002); we attribute the discrepancy 
to the different modeling of the continuum emission (Read & Stevens 2002 did not include 
pn data, neither considered the contribution from point sources). However, it is worth of 
mentioning that both sets of abundances indicate a low O abundance compared to the other 
a-elements. 

To further check the consistency and robustness of our estimates, we repeated the fit two 
times considering only the RGS data and using i) the above described thermal plus power-law 
model and ii) a, thermal model alone. This last check, in particular, most closely matches the 
analysis previously made by Read & Stevens (2002) to the point that the discrepancies are 
within the errors and can be attributed to the differences in the used instrument calibration. 
The abundances derived from the three different models are plotted in Fig. 5. Those obtained 
combining pn and RGS have significantly smaller uncertainties. One can sec that different 
models have mainly the effect to almost rigidly scale all the abundances, with minor impact 
on the abundance ratios. 

The low Oxygen abundance is somehow difficult to explain in the framework of a- 

element enhancement by type II SN explosions. Thus, to check its reliability we included 
in the fit also the channels corresponding to the 20-23 A interval, where the OVII triplet 
is observed (although with very few counts). We find no significant difference neither in 
Oxygen abundance nor in DEM distribution (the OVIII A19/0VII A22 ratio is extremely 
sensitive to the plasma temperature). Ne is also somewhat under- abundant. We checked for 
possible instrumental effects, such as variations in the RGS effective area in proximity of the 
Ne and O lines; however the effective area in this region is a rather smooth function of the 
wavelength, so that there should be no instrumental issues and these low abundances are 
likely to be real. 

Since M82 is observed edge — on, the X-ray emission of its nuclear region is subject to 
quite heavy absorption which could affect the O lines detected at low energies. However, 
as shown in Fig. 6, the contour plot of the Oxygen abundance towards the foreground 
absorption, suggest a rather weak correlation. In the next section we discuss other physical 
possibilities to explain the O under-abundance issue. 
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Fig. 4. — 68. 90 and 99% confidence contours of the Iron abundance vs intensity of the 
thermal component. For easy of reading the y-axis is shown both in hnear and logarithmic 
scale. 



- 15 - 



5 



o 
X 



X 



0,5 



0.1 




0.6 



0.4 



0.2 







-0.2 



-0.4 



-0.6 



-0.8 



X 

X 



Fe 







Ne 



Mg 



Si 



Fig. 5. — Best-fit clement abundances of the hot gas for different models. Filled circles: fit 
of both RGS and pn data, thermal plus power-law model. Open squares: RGS data alone, 
thermal plus power-law model. Filled triangles: RGS data alone, thermal model alone. For 
easy of reading the y-axis is shown both in linear and logarithmic scale. 
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Fig. 6. — Hot gas-phase Oxygen abundance vs absorption. For easy of reading the y-axis is 
shown both in hnear and logarithmic scale. The correlation, if any, is rather weak, and the 
Oxygen abundance can be well constrained. 
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Fig. 7. — Stellar, hot and cold gas abundance ratios relative to Iron. Cold gas-phase abun- 
dances are from different data sets in the literature, hence they are not homogeneously 
determined. The dashed hne indicates the solar values. For easy of reading the y-axis is 
shown both in linear and logarithmic scale. 
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5. Discussion 

In M82 all the three components, namely hot and cold gas-phases and stars trace a very 
similar Iron abundance, the average value being [Fe/H]~ —0.35 dex. Indeed, since Iron is 
mainly produced by SNe la, it is expected to be released in the ISM only after ~ 1 Gyr from 
the local onset of SF. 

At variance, a— elements (O, Ne, Mg, Si, S, Ca, Ti) are predominantly released by 
SNe II with massive progenitors on much shorter timescales. 

Stars trace an average [<Si,Mg,Ca>/Fe]~ 0.4 dex, while the hot gas suggests an average 
[<Si,Mg,S>/Fe]~ 0.5 dex. Such an overall [a/Fe] enhancement in M82 is fully consistent 
with a standard chemical evolution scenario, where the ISM in the nuclear regions of massive 
SB galaxies is mainly enriched by the products of SNe II explosions (see e.g. Arnett 1995 for 
a review) occurring in recursive bursts of SF of relatively short duration. 

Fig. 7 shows the various abundance ratios relative to Iron, as measured in stars, hot 
and cold gas-phases. Since the cold gas-phase abundances are from different data sets in 
the literature, hence not homogeneously determined, the inferred abundance ratios should 
be regarded with caution. 

When comparing the hot and cold gas-phase abundances, one finds Ne over-abundant. 
Si and especially S significantly under-abundant in the cold gas. Forster Schreiber et al. 
(2001) suggests S and to a lower level Si depletion onto interstellar diist grains. However, 
other metals like Fe in particular, with even higher degree of incorporation into dust (Savage 
& Sembach 1996), should be severely depleted in the cold gas, which does not seem the case. 
A possible explanation is that the optical lines used to infer Fe and O abundances (Alloin, 
Collin-Souffrin, & Joly 1979) and the the mid/far IR lines used to infer Ne and S abundances 
(Lord et al. 2001; Forster Schreiber et al. 2001) trace different nebular sub-structures within 
the central few hundreds pc, somewhat chemically dishomogeneous or with different dust 
content. 

The [0/Fe] abundance ratio in M82 is even more puzzling to interpret. Before the onset 
of the current burst of SF, the ISM was enriched in O as well as in the other a-elements, as 
suggested by our stellar and cold gas-phase abundances. The O under- abundance measured 
in the hot gas cannot be modeled with a standard nucleosynthesis from the present gen- 
eration of SNe II. Umeda et al. (2002) suggested explosive nucleosynthesis in core-collapse 
hypernovae but they predicted a large under- abundance of Ne and Mg as well, which is not 
observed. [0/Fe] under-abundance can be also explained with a major Fe enrichment by 
current SN-Ia explosions from previous generations of stars, but in this case all the [a/Fe] 
ratios should be low, again contrary to what observed. O might be also locked into dust 
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grains. Indeed, SNe II are known to produce significant amounts of dust (Todini & Ferrara 
2001) and O, being a basic constituent of dust grains, is expected to suffer of severe dust 
depletion. However, otlier metals like Fc in particular, but also Si and Mg should be de- 
pleted into dust (Savage & Sembach 1996), leaving almost unchanged or eventually further 
enhancing the [a/Fe] ratio. 

6. Conclusions 

Our abundance analysis of the stellar and hot gas-phase components in the nuclear 
region of M82 indicate an Iron abundance about half Solar and an overall a— enhancement 
by a factor between 2 and 3 with respect to Iron. 

These abundance patterns can be easily explained within a standard nucleosynthesis 
scenario where the ISM is mainly enriched by SNe II on relative short timescales and with 
a star formation process occurring in recursive bursts. 

Oxygen behaves in a strange fashion. It is over-abundant in stars and cold gas, similarly 
to the other a-elements, while is significantly under- abundant in the hot gas. Major calibra- 
tions and/or modeling problems seem unlikely. Hypernovae nucleosynthesis, dust depletion, 
SN-Ia enrichment can somehow explain an Oxygen under-abundance but other metals should 
follow the same pattern, while they do not. 

Somewhat exotic threshold effects and/or depletion mechanisms preferentially affecting 
Oxygen could be at work in the nuclear region of M82, but presently this issue remains 
controversial. 
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